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^ ABSTRACT 

o 

We present radio observations of a region in the vicinity of the young pulsar 
PSRB1853+01 in the supernova remnant W44. The pulsar is located at the apex 
Qh of an extended feature with cometary morphology. We argue on the basis of its 

Q morphology and its spectral index and polarization properties that this is a synchrotron 

nebula produced by the spin down energy of the pulsar. The geometry and physical 
Ci parameters of this pulsar-powered nebula and W 44 are used to derive three different 

measures of the pulsar's transverse velocity. A range of estimates between 315 and 
470 km s~^ are derived, resulting in a typical value of 375 km s~^. The observed 
synchrotron spectrum from radio to X-ray wavelengths is used to put constraints on 
the energetics of the nebula and to derive the parameters of the pulsar wind. 



Subject headings: Pulsars: individual (PSR B1853-|-01) - supernova remnants — ISM: 
individual (W44) 
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1. Introduction 

W44 (G 34.7—0.4) belongs to a small group of supernova remnants in our Galaxy which 
contain a young, active pulsar. At radio wavelengths W44 has a shell-type morphology, about 30' 
in diameter, with the bulk of its radio emission concentrated in a series of knots and filaments 
(Jones, Smith & Angellini 1993). Its X-ray emission, mostly thermal in origin, is confined 
within the radio shell and is centrally-peaked (Smith et al. 1985, Rho et al. 1994). Although 
PSR B1853-I-01 is offset nearly 9' south of the center of W44, the similarities in their age and 
distance make a compelling case for a real, physical association between them (Wolszczan, Cordes 
k Dewey 1991). PSRB1853-F01 has a period P of 267 ms and a period derivative P of 208x10"^^ s 
s~^, implying a surface dipole field Bq of 7.6x10^^ G and rotational energy loss rate E of 4.3x10"^^ 
erg s~^ (Wolszczan et al. 1991). Its small characteristic spindown age Tc of only 20,000 yrs makes 
PSR B1853-I-01 one of the 10 youngest pulsars in our Galaxy (Frail, Goss & Whiteoak 1994). 

It is widely thought that young pulsars efficiently transfer their rotational spindown energy E 
into a relativistic wind with some mean Lorentz factor 7 (Michel 1969, Rees & Gunn 1974, Kennel 
et al. 1983). The energy in the wind is divided between magnetic and electric fiux (i.e. Poynting 
fiux) and particle energy fiux, parametrized by the ratio a. From current models of pulsar 
magnetospheres it is expected that the wind will be magnetically dominated {a ^1) out to at least 
the light cylinder radius (Arons 1992), but self-consistent models of the Crab nebula and pulsar 
require that a further out (Kennel & Coroniti 1984, Emmering & Chevalier 1987). Several 
solutions have been proposed to explain this discrepancy (Begelman & Li 1994, Michel 1994, 
Melatos & Melrose 1996), including that by Coroniti (1990), who invoked magnetic reconnection 
from field reversals to lower a by taking energy from the B-field to accelerate particles. It is not 
known how common this behavior might be, since apart from the Crab pulsar and PSR B1957-|-20 
(Kulkarni et al. 1992), there are very few measurements of ci. Another wind parameter of interest 
is the particle composition, expressed as k where k=U8/Ue, a ratio between the energy density 
in electrons Ue versus all other particles Ui (including positrons, protons and heavy ions). While 
most models assume an electron-positron pair plasma (k=l), an electron-ion plasma (k^l) is not 
ruled out, and may in fact be required to produce the observed power law distribution of the 
post-shock particles (Hoshino et al. 1992, Gallant et al. 1992). 

The best way to constrain the wind parameters is through the study of synchrotron nebulae 
which appear as diffuse sources in the vicinity of some pulsars. The Crab nebula and other 
"plerions" are the best-known examples (Weiler & Sramek 1988). However, there are radio "wisps" 
and similar elongated structures, including the radio emission near the Vela pulsar (Bietenholz, 
Frail k Hankins 1991), the Crab pulsar (Bietenholz k Kronberg 1992), the CTB 80 pulsar 
(Strom 1987), the X-ray point source in 3C58 (Frail k Moffett 1993), and the source surrounding 
PSRB1757— 24 (Frail k Kulkarni 1991). These features are thought to originate at the termination 
shock from a pulsar wind, formed where the wind first interacts with the surrounding medium. 
Thus the investigation of pulsar wind nebulae, their morphology, polarization structure, and 
spectral and temporal behavior, can put constraints on existing theoretical models. 
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Recently, Jones et al. (1993) reported a detailed X-ray and radio study of W44. In their 
20-cm VLA radio images, they claim the detection of two discrete sources in the vicinity of the 
pulsar: weak extended radio emission at 20 cm, 12" north of the pulsar with a flux of 2 mJy, and 
a "diamond-shaped" nebula 1'.5 north, with a flux of 60 mJy. Jones et al. (1993) proposed that 
this diamond-shaped radio feature was a pulsar-powered synchrotron nebula. However, doubts 
remained about this identiflcation because of the large offset between the pulsar and the nebula 
and the uncertain nature of the extended emission 12" away from the pulsar. The present study 
was undertaken to address these difficulties by providing an improved position for PSR B1853-|-01 
and imaging the region around the pulsar at a variety of radio frequencies. 



2. Observations and Results 

The interferometric position of the pulsar was determined from two datasets; the flrst was 
archival data taken on 1990 June 18 with the VLA in its BnA array conflguration, and the second 
was on 1994 June 12 in the B array. Both observations used the right circular hands of polarization 
at 1465 MHz and 1515 MHz. In B array data there is only one radio source within a radius of 45" 
of the nominal timing position from Wolszczan et al. (1991). A flt to the peak gives Q!(1950) = 18'' 
53™ 38. 33^ ^(1950) = -F01° 09' 26.1". This position is approximately 4" to the northwest of the 
timing position. Since the timing position is based on only 18 months of timing data, and the 
pulsar exhibits timing noise comparable to that of the Crab pulsar (Wolszczan et al. 1991) such a 
discrepancy is not unexpected. At the higher resolution of the BnA array data (beamwidth of 1"6 
versus 4") the emission in this region appears resolved, with three equally bright compact features 
clustered at the location of the B array peak. It appears that during the BnA array observations 
PSR B1853-I-01 was anomalously weak (<0.3 mJy) and thus we cannot uniquely identify which of 
these peaks is the pulsar and which are partially resolved extended emission. While the formal 
errors on the B array flts are of order ±0.15", we assign an error of ±1.5" to reflect the uncertainty 
in correctly identifying the pulsar in the presence of extended emission. 

The extended emission in the vicinity of the pulsar was imaged at 4860 MHz and 8440 MHz 
on 1995 May 20 with the VLA in its D array conflguration. Two 50 MHz channels were measured 
at 6 cm and 3.6 cm in both parallel and cross hands of circular polarization. All data reduction 
and calibration was done following standard practice in use at the VLA. In Figure lb we show the 
8.4 GHz (A=3.6 cm) image of the pulsar, and for perspective in Figure la we include a 1.4 GHz 
image showing the location of the pulsar in relation to the larger W44 remnant. The 1.4 GHz 
image is part of a comprehensive radio. X-ray and optical study of the W44 remnant (Giacani et 
al. 1996). Unlike a similar image at 1.4 GHz by Jones et al. (1993), care was taken to add the 
missing short order spacings using single dish data. Consequently the true nature of the emission 
around PSRB1853±01 is revealed. 

With the addition of short spacings the two discrete sources seen by Jones et al. (1993), the 
"diamond" nebula and the extended emission 12" north of the pulsar, form a single structure. 
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The morphology is even better defined in Fig. lb. The pulsar is located at the apex of a 
cometary-shaped radio source 2'.5 in extent. A narrow "contrail" leads away from the pulsar 
to the northwest where it widens out and the continuum emission peaks. There is a "kink" in 
the nebula about 30" from PSRB1853+01, that could be produced by density or temperature 
variations within W44, and should be detectable as fiuctuations in the X-ray emissivity. Similar 
axisymmetric structures have been seen at radio wavelengths including G 5.4—1.2 (Frail & Kulkarni 
1991), G 359.2-0.8 (Yusef-Zadeh & Bally 1987, Predehl & Kulkarni 1994) and G 354.1+0.1 (Frail 
et al. 1994). Like the nebula in Fig. 1, these are best understood as synchrotron nebulae that 
arise as a result of the spin down of young, energetic pulsars. 

Several other lines of evidence support the contention, first put forth by Jones et al. (1993), 
that this is a synchrotron nebula generated by the pulsar and not just a superposition of a filament 
from W44 projected onto PSRB1853+01. The total fiux density of the nebular emission around 
PSRB1853+01 is 287±126 mJy at 20 cm and 220±26 mJy at 3.6 cm. A radio spectrum between 
330 MHz and 8.4 GHz using data from Giacani et al. (1996) is given in Figure 2. The larger error 
bars for the lower frequency radio data refiect the difficulty in separating this emission from its 
surroundings. At these frequencies this nebular emission is < 10""^ of the total fiux density of 
W44 (Kovalenko, Pynzar' & Udal'tsov 1994a, b). We include an X-ray upper limit from the Ginga 
satellite determined by Harrus, Hughes & Helfand (1996) who also report on the detection of an 
extended, hard X-ray source near this fiux level with the ASCA satellite. A weighted least-squares 
fit to the radio portion of the spectrum yields cir = —0.12 ± 0.04 (where Sj, oc /^"^), substantially 
fiatter than cir = —0.33 for the rest of W44 (Kovalenko et al. 1994a, b). The emission is polarized 
with a degree of polarization at 3.6 cm and 6 cm of approximately 17±4%. Although W44 exhibits 
a high degree of polarization of up to 20% at frequencies above 10 GHz (Kundu & Velusamy 
1972), the emission is strongly depolarized at lower frequencies, dropping to a peak of 3% at 2.7 
GHz (Velusamy & Kundu 1974). At 5 GHz (6 cm) we estimate from the data of Whiteoak & 
Gardner (1971) that the degree of polarization for W44 in the vicinity of the pulsar is only 3.5%. 
At this same frequency our data shows that the polarized fiux from the compact nebula stands 
out from the body of the surrounding emission, confirming its much higher degree of polarization 
compared to W 44 as a whole. 

Thus the case for the synchrotron source around PSR B1853-|-01 being a pulsar-powered wind 
nebula (PWN) rests on its unusual cometary morphology with the pulsar at its apex, and the 
fact that its radio spectrum and degree of polarization distinguish it from the rest of W44. In 
these latter two properties this PWN has much in common with other pulsar-powered nebula or 
"plerions" of which the Crab nebula is the prototype (Weiler & Sramek 1988). In the next section 
we use the observed properties of the PWN to put constraints on the pulsar parameters. 
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The morphology of the PWN can be used to put independent constraints on both the 
magnitude and direction of the transverse velocity VpsR of PSRB1853+01. Wolszczan et al. 
(1991) and Frail et al. (1994) earlier derived Vpsr of PSRB1853+01 on the basis of the pulsar's 
displacement from the geometric center of the remnant and its age. A better estimate of these can 
be made, since the birthplace of PSRB1853+01 can be located by three different methods. The 
geometric center of the radio continuum shell (Jones et al. 1993), the dynamical center of the HI 
shell (Koo & Heiles 1995), and the peak of the X-ray emission (ROSAT: Rho et al. 1994, Einstein: 
Seward 1990) all pinpoint the probable center and agree quite well, giving an angular displacement 
of PSR 61853+01 from this point of 8f6±0f5. This point lies along a line that bisects the axis of 
symmetry of the PWN, suggesting that the pulsar originated from this location. It follows that 
the implied velocity of PSR B1853+01 is VpsR = 370^3^20"^ km s~^ (±25 km s~^), where the age 
for W44 of 20,000^20 years is from Tc of PSRB1853+01, which is likely an upper limit. 

Another estimate of the velocity of PSRB1853+01 is possible. For any supernova remnant 
we can express the shock velocity at a time t as = Cq Rs/^, where R^ is the radius and Cq is a 
constant equal to 2/5 for a remnant in the Sedov phase (ShuU et al. 1989). Similarly, the pulsar's 
velocity is VpsR = /3Rs/t, where (3 is the fractional displacement of the pulsar from the center of 
the remnant. Combining these two equations we see that the pulsar velocity and shock velocity 
are related by Vpsr=Vs /3/co. For PSRB1853+01 we estimate /3 ~0.5, while Koo & Heiles (1995) 
derive Vs=330 km s~^ for the radio continuum shell of W44. If W44 is near the Sedov phase (as 
favored by Koo & Heiles) then for Cq — 0.35-0.45 we derive VpsR=370-470 km s~^. 

A third way to estimate the pulsar velocity is through the pressure balance condition that 
must exist between the PWN and the hot gas interior to W44. According to Rees & Gunn (1974) 
the relativistic pulsar wind flows outward until it terminates at a point r^, where it comes into 
pressure equilibrium with the surrounding medium. For W44 this pressure is given by 

— = 5.7 X 10-^°6'2d3-2 erg cm'^ 

where we have assumed an isotropic wind and 9^ is the angular radius (in arcsec) where the 
reverse shock is formed. This region appears unresolved in Fig. 1 but from our high resolution 
data we estimate that 9^ <1". The distance of Sd^ kpc is well-established by HI absorption and 
1720 MHz maser lines (Caswell et al. 1975, Claussen et al. 1996). Some fraction of E could be 
converted into high-energy radiation rather than particles and held energy in the wind but upper 
limits for the detection of pulsed emission from PSR B1853-|-01 with ROSAT (0.1-2.4 keV, Rho et 
al. 1994) and EGRET (>100 Mev, Thompson et al. 1994a) suggest that the production efficiency 
T] is low (<10%). The unusual cometary morphology leaves little doubt that the source of external 
pressure conflning the flow is the ram pressure pVp^^ due to the pulsar's motion through W44, 

pY^sR = 2.3 X lO-^NoV^oo erg cm'^. (2) 



(1) 
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We have allowed for a 10% composition of He by mass in calculating the density no. Vioo is the 
pulsar velocity in units of 100 km s~^. Equating eqns (1) and (2) and using no=0.25 cm~"^ near 
PSRB1853+01 from Jones et al. (1993) we derive Vpsr > 3156'jM3-^ km s'^ 

Beyond r^ the wind becomes sub-relativistic and is swept back and kept separate from the 
shocked external gas by a bow shock that leads ahead of the pulsar (Cheng 1983, Wang, Li & 
Begelman 1993). A synchrotron trail forms in a direction opposite to the pulsar's motion, confined 
by the external medium with a length determined by adiabatic and synchrotron losses. For the 2'.5 
long PWN trailing behind PSRB1853+01 the thermal pressure from the hot X-ray emitting gas 
in W44 is approximately 6-8x10"^° erg cm~"^ (Jones et al. 1993), comparable to the ram pressure 
in eqn. 2 and several orders of magnitude above the pressure in the general interstellar medium 
of 4xlO~^^(no To/3000) erg cm~^ (Kulkarni & Heiles 1988). This trail points in a northwest 
direction, supporting our contention made earlier, that the pulsar originated close to the geometric 
center of W44. 

In summary, we find that all three independent methods yield an estimate for VpsR which 
agree remarkably well, suggesting that in the absence of proper motion measurements, inferred 
velocities are basically sound. We adopt a value VpsR=375 km s~^ for the transverse velocity of 
PSR B1853-I-01 and predict proper motion of order 25 mas yr~^. 

The observed synchrotron spectrum from radio to X-ray wavelengths (Fig. 2) can be used 
to put constraints on the energetics of the nebula (Helfand & Becker 1987, Salter et al. 1989) 
and the parameters of the pulsar wind discussed in §1 (i.e. 7, a and k). From the spectrum 
displayed in Fig. 2, it can be seen that the Ginga X-ray upper limit does not lie on a smooth 
extrapolation of the radio spectrum, implying that the spectrum must break at some z/p. Spectral 
breaks hold the key to understanding the physical conditions of the nebula and the pulsar wind. 
If the break results from synchrotron losses then, given the age of the source t^es, we can estimate 
the nebular field B„. From VpsR and the linear size of the PWN, we derive t^es — 5700 yrs 
and from the standard formulae in Pacholczyk (1970) B„ = 1040//G (//p/lO^^ Hz)-i/3(tres/1000 
yrs)"^/"^. A Crab-like spectral index for the X-ray emission ax = — 1 gives z/p = 5.7 X 10^^ Hz and 
Bn — 400 //G. If ax is steep, say —1.5, then z/p = 1.2 X 10^"* Hz and B„ ~70//G. An upper limit 
of ax = —0.81 is derived by setting z/p to 10^° Hz, corresponding to the far end of the observed 
radio spectrum. In this case B„ ~1500//G. 

A lower limit to the value of z/p can be derived by using other measures of the magnetic 
field. The maximum field reached in the pulsar wind is when it remains Poynting fiux-dominated 
(a > 1) out to the reverse shock. Then from eqn. 1 we can estimate the equivalent magnetic 
field Bp = ^2E/rg c which would contain the energy density of the wind at r^ (Thompson et al. 
1994b). Similarly, if we assume that the magnetic field, which is Bo at the surface of the neutron 
star at a radius R, falls off as a dipole out to the light cylinder radius r/c = il/c and is thereafter 
toroidal to r^, then Bs=Bo(R/r;c)"^(r;c/i's). Both values give B^ ~110 //G but it could be much 
lower if a decreases from r/c to r^, as seems necessary for the Crab nebula (Kennel & Coroniti 
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1984, Emmering & Chevalier 1987). Some enhancement of the postshock B„ is expected over 
but at least in the case of a strong relativistic shock, B„ < SB^ (Kennel & Coroniti 1984). From 
the values of B„ derived earlier this allows us to limit z/b > 10^^ Hz. 

It follows that the Lorentz factor 7 for the electrons responsible for the emission near z/b is 
7 ~ 10^(//b/10^2h2)1/2 (B/100//G)-^/2 (Pacholczyk 1970). Larger values of 7 and its evolution 
away from the pulsar are measurable in principle from high resolution X-ray observations. 
Integrating the spectrum from 10'' Hz to z/b = 10^^ Hz we derive a radio luminosity Lr of 
1.5xlO"^"^d3 ergs s~^. Likewise from the Ginga data the X-ray luminosity Lx is 1.3xlO"^"^d3 ergs 
s~^. Both Lr and Lx are only ~0.003E. The empirical relation of Seward & Wang (1988) between 
Lx and E for PWN predicts Lx ~0.002E for the nebula around PSR B1853-F01. Adopting the 
equipartition condition (Pacholczyk 1970) we derive Beg > 60 //G (l-|-k)^/'' (z/r > 10^^ Hz) and a 
minimum energy E^m > 1.3 X 10"*^ ergs (l-|-k)'*/'' for the PWN. It is not hard to show that for a 
dipole rotator the energy deposited by the pulsar into the nebula is Ext^es [t2o/(t20— Ves)], and 
thus an upper limit on k can be derived by requiring E^m — 1.4Extres— 10'*'' ergs. At z/b = 10*^ 
Hz k=30 but it falls to k=5 at z/b = lO*"* Hz, suggesting that, while ions may exist in the wind, 
they are evidently energetically less important than some models assume (Gallant & Arons 1994). 

Finally we comment that if the magnetic field B„ were better known (from z/b) then the 
energy density ratio between magnetic fields Ub and particles Up in the nebula =Ub/Up could 

7 /2 

be determined (since (T„ ocVB„ /(l-|-k)LR, where V is the nebular volume), leading ultimately to 
a of the wind. The constraints on the pulsar wind diagnostics deduced here will improve once a 
better measurement of z/r is made and the data are interpreted with more complex models (e.g. 
Reynolds & Chevalier 1984, Reynolds & Chanan 1984) which will need to incorporate continuous 
particle/field injection into the cometary-shaped nebula produced by the pulsar's space motion. 

The Very Large Array (VLA) is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. This research has made use of the Simbad 
database, operated at CDS, Strasbourg, France. DAE thanks S. Kulkarni, J. Navarro and G. 
Vasisht for useful discussions. 
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Fig. 1. — (a) A 20-cm radio continuum image of the supernova remnant W44. The inset region 
is shown in more detail in (b) with a 3.6 cm image showing the pulsar wind nebula around 
PSRB1853+01. The position of the pulsar is indicated by a cross ( + ). Contour intervals are —0.6, 
0.6, 1, 1.4, 1.8, 2.2, 2.6, 3, 4 mJy beam~^. The rms noise is 0.2 mJy beam~^. The synthesized 
beam size is 8'.'9x7'.'8, P.A.=-24°. 

Fig. 2. — The spectrum of the pulsar wind nebula around PSRB1853+01 over the range from 10® 
Hz to 10^® Hz. Solid line with slope cir in the radio regime is a weighted least-squares fit to the 
data points discussed in the text. The X-ray point is an upper limit from Ginga data (see text). A 
Crab-like power-law spectral index ax = —1.5 is drawn through the X-ray data point and intersects 
the extrapolation of the radio spectrum at a break frequency z/b . 



